Common recognition topology of mex transporters of Pseudomonas aeruginosa revealed by molecular modelling

The secondary transporters of the resistance-nodulation-cell division (RND) superfamily mediate multidrug resistance in Gram-negative bacteria like Pseudomonas aeruginosa. Among these RND transporters, MexB, MexF, and MexY, with partly overlapping specificities, have been implicated in pathogenicity. Only the structure of the former has been resolved experimentally, which together with the lack of data about the functional dynamics of the full set of transporters, limited a systematic investigation of the molecular determinants defining their peculiar and shared features. In a previous work (Ramaswamy et al., Front. Microbiol., 2018, 9, 1144), we compared at an atomistic level the two main putative recognition sites (named access and deep binding pockets) of MexB and MexY. In this work, we expand the comparison by performing extended molecular dynamics (MD) simulations of these transporters and the pathologically relevant transporter MexF. We employed a more realistic model of the inner phospholipid membrane of P. aeruginosa and more accurate force-fields. To elucidate structure/dynamics-activity relationships we performed physico-chemical analyses and mapped the binding propensities of several organic probes on all transporters. Our data revealed the presence, also in MexF, of a few multifunctional sites at locations equivalent to the access and deep binding pockets detected in MexB. Furthermore, we report for the first time about the multidrug binding abilities of two out of five gates of the channels deputed to peripheral (early) recognition of substrates. Overall, our findings help to define a common “recognition topology” characterizing Mex transporters, which can be exploited to optimize transport and inhibition propensities of antimicrobial compounds

Tripartite efflux pumps of the RND superfamily: what did we learn from computational studies?

Bacterial resistance to antibiotics has been long recognized as a priority to address for human health. Among all micro-organisms, the so-called multi -drug resistant (MDR) bacteria, which are resistant to most, if not all drugs in our current arsenal, are particularly worrisome. The World Health Organization has prioritized the ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Enterobacter species) pathogens, which include four Gram-negative bacterial species. In these bacteria, active extrusion of antimicrobial compounds out of the cell by means of 'molecular guns' known as efflux pumps is a main determinant of MDR phenotypes. The resistance-nodulation- cell division (RND) superfamily of efflux pumps connecting the inner and outer membrane in Gram-negative bacteria is crucial to the onset of MDR and virulence, as well as biofilm formation. Thus, understanding the molecular basis of the interaction of antibiotics and inhibitors with these pumps is key to the design of more effective therapeutics. With the aim to contribute to this challenge, and complement and inspire experimental research, in silico studies on RND efflux pumps have flourished in recent decades. Here, we review a selection of such investigations addressing the main determinants behind the polyspecificity of these pumps, the mechanisms of substrate recognition, transport and inhibition, as well as the relevance of their assembly for proper functioning, and the role of protein-lipid interactions. The journey will end with a perspective on the role of computer simulations in addressing the challenges posed by these beautifully complex machineries and in supporting the fight against the spread of MDR bacteria

Molecular electrometer and binding of cations to phospholipid bilayers

Despite the vast amount of experimental and theoretical studies on the binding affinity of cations -especially the biologically relevant Na+ and Ca2+ - for phospholipid bilayers, there is no consensus in the literature. Here we show that by interpreting changes in the choline headgroup order parameters according to the 'molecular electrometer' concept [Seelig et al., Biochemistry, 1987, 26, 7535], one can directly compare the ion binding affinities between simulations and experiments. Our findings strongly support the view that in contrast to Ca2+ and other multivalent ions, Na+ and other monovalent ions (except Li+) do not specifically bind to phosphatidylcholine lipid bilayers at sub-molar concentrations. However, the Na+ binding affinity was overestimated by several molecular dynamics simulation models, resulting in artificially positively charged bilayers and exaggerated structural effects in the lipid headgroups. While qualitatively correct headgroup order parameter response was observed with Ca2+ binding in all the tested models, no model had sufficient quantitative accuracy to interpret the Ca2+: lipid stoichiometry or the induced atomistic resolution structural changes. All scientific contributions to this open collaboration work were made publicly, using nmrlipids. blogspot.fi as the main communication platform.Peer reviewe

Novel Changes in Discoidal High Density Lipoprotein Morphology: A Molecular Dynamics Study

AbstractApoA-I is a uniquely flexible lipid-scavenging protein capable of incorporating phospholipids into stable particles. Here we report molecular dynamics simulations on a series of progressively smaller discoidal high density lipoprotein particles produced by incremental removal of palmitoyloleoylphosphatidylcholine via four different pathways. The starting model contained 160 palmitoyloleoylphosphatidylcholines and a belt of two antiparallel amphipathic helical lipid-associating domains of apolipoprotein (apo) A-I. The results are particularly compelling. After a few nanoseconds of molecular dynamics simulation, independent of the starting particle and method of size reduction, all simulated double belts of the four lipidated apoA-I particles have helical domains that impressively approximate the x-ray crystal structure of lipid-free apoA-I, particularly between residues 88 and 186. These results provide atomic resolution models for two of the particles produced by in vitro reconstitution of nascent high density lipoprotein particles. These particles, measuring 95Å and 78Å by nondenaturing gradient gel electrophoresis, correspond in composition and in size/shape (by negative stain electron microscopy) to the simulated particles with molar ratios of 100:2 and 50:2, respectively. The lipids of the 100:2 particle family form minimal surfaces at their monolayer-monolayer interface, whereas the 50:2 particle family displays a lipid pocket capable of binding a dynamic range of phospholipid molecules

Role of Lipids in Spheroidal High Density Lipoproteins

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly

Partial density profiles of lipids and ions

<p>Partial density profiles of lipids and ions calculated using the script available on Zenodo. Each simulated system contains 600 DPPC molecules, 150 mM NaCl and 30 waters/lipid. The ion model is the one described by Roux. The analysis has been performed over the last 40 ns of each 100 ns trajectory. All simulations have been performed with the Slipids force field.</p

Application of molecular modelling and EPR spectroscopy to lipid membranes - a combined approach

Knowledge of molecular interactions, thermodynamics, temperature and system composition effects are crucial for understanding the role that different lipids play in vital life processes in biological membranes. This knowledge is also important for understanding the impact that electromagnetic fields have on the order and mobility of molecules in lipid bilayers. The last decade has seen radical improvement in the molecular modelling of complex molecular and bio-molecular systems including lipid bilayers using Molecular Dynamics (MD) simulation techniques. MD simulations are now much faster and more accurate allowing researchers to predict complex molecularphenomena using actual structures.In this paper we present our recent results on the application of large scale MD simulations to phospholipid bilayers under different composition and conditions.Examples include: both all-atom and coarse-grain large scale MD simulations of binary and ternary compositions of lipid bilayers, modelling of the penetration of gas molecules (O2) in lipid bilayers, the effects of antimicrobial peptides on biological membranes, separation of lipid microdomains as a model for the study of lipid rafts. We also report MD simulations on lipid bilayers doped with structurally different nitroxide spin probes that are employed inexperimental variable temperature EPR spectroscopy. Finally, our recent preliminary results of all-atom MD modelling of the lipid bilayers subjected to microwave electric fields are also presented and discussed

Dynamics of Activation of Lecithin: Cholesterol Acyltransferase by Apolipoprotein A-I

The product of transesterification of phospholipid acyl chains and unesterified cholesterol (UC) by the enzyme lecithin:cholesterol acyltransferase (LCAT) is cholesteryl ester (CE). Activation of LCAT by apolipoprotein (apo) A-I on nascent (discoidal) high-density lipoproteins (HDL) is essential for formation of mature (spheroidal) HDL during the antiatherogenic process of reverse cholesterol transport. Here we report all-atom and coarse-grained (CG) molecular dynamics (MD) simulations of HDL particles that have major implications for mechanisms of I-CAT activation. Both the all-atom and CG simulations provide support for a model in which the helix 5/5 domains or apoA-I create an amphipathic "presentation tunnel" that exposes methyl ends of acyl chains at the bilayer center to solvent. Further, CG simulations show that UC also becomes inserted with high efficiency into the amphipathic presentation tunnel with its hydroxyl moiety (UC-OH) exposed to solvent; these results are consistent with trajectory analyses of the all-atom simulations showing that UC is being concentrated in the vicinity of the presentation tunnel. Finally, consistent with known product inhibition of CE-rich HDL by CE, CG simulations of CE-rich spheroidal HDL indicate partial blockage of the amphipathic presentation tunnel by CE. These results lead us to propose the following working hypothesis. After attachment of LCAT to discoidal HDL, the helix 5/5 domains in apoA-I form amphipathic presentation tunnels for migration of hydrophobic acyl chains and amphipathic UC from the bilayer to the phospholipase A2-like and esterification active sites of LCAT, respectively. This hypothesis is currently being tested by site-directed mutagenesis

Molecular dynamics simulations of monomeric apolipoprotein A-I from a recent X-ray structure

We have examined the X-ray crystal structure recently refined by Ajees and colleagues (Ajees et al. 2006) for monomeric apolipoprotein A-I (apoA-I). Because the structure, which has been crystallized together with chromium organic compounds, possesses a substantially higher percentage of alpha helicity than is generally estimated experimentally for the lipid-free monomeric apoA-I in solution (~80% vs ~50%), we have performed molecular dynamics (MD) simulations for ~10 ns of the model in order to explore the dynamic behavior of the single apoA-I monomer at a physiological salt concentration and a temperature range of 310-410 K. While 10 ns simulation is only a starting point, a few important observations have been made: i) the percentage of alpha helicity decreased substantially to below 70% (i.e., towards a lower experimental estimate); ii) the structure became more globular in overall appearance; iii) the flexible N-terminal domain (amino acid residues 1 to 43) has lost most of its alpha helicity; iv) the hydrophobic core of the 4-helix bundle is defined by stacking of a cluster of aromatic amino acid residues, outlined by a shell of aliphatic hydrophobic residues; v) The four helix bundle portion of the simulated structure is clustering around a pronounced stacking of aromatic residues derived from all four helixes and the aromatic cluster is overlaid by a shell of aliphatic hydrophobic residues. We conjecture that this aromatic cluster and its surrounding hydrophobic residues are the driving force for creation of a dynamic (molten globular) four helix bundle arrangement in lipid-free monomeric apoA-I in solution. This work was supported by NIH grant